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ABSTRACT 

Context. The Cygnus region harbours a huge complex of massive stars at a distance of 1.0-2.0kpc from us. About 170 O stars are 
distributed over several OB associations, among which the Cyg OB2 cluster is by far the largest with about 100-120 O stars. As a 
consequence of their successive nuclear-burning episodes, these massive stars inject large quantities of radioactive nuclei into the 
interstellar medium such as 26 Al and 60 Fe. The gamma-ray line signal from the latter is a solid tracer of ongoing nucleosynthesis. 
Aims. We want to compare the decay emission from the Cygnus region with predictions of recently improved stellar models. As a first 
step, we establish observational constraints upon the gamma-ray line emission from 26 Al and 60 Fe, with particular emphasis placed 
on separating the emission due to the Cygnus complex from the foreground and background mean Galactic contributions. 
Methods. We used the high-resolution gamma-ray spectrometer INTEGRAL/SPI to analyse the 26 Al and 60 Fe decay signal from the 
Cygnus region. The weak gamma-ray line emissions at 1809 and 1 173/1332 keV have been characterised in terms of photometry, 
spectrometry, and source morphology. 

Results. The 1809 keV emission from Cygnus is centred on the position of the Cyg OB2 cluster and has a typical size of 9° or more. 
The total 1809keV flux from the Cygnus region is (6.0 ± 1.0) x 10~ 5 phcnT 2 s , but the flux really attributable to the Cygnus complex 
reduces to (3.9± l.l)x 10 s phcnT 2 s . The 1809 keV line centroid is agrees with expectations considering the direction and distance 
of the Cygnus complex, and the line width is consistent with typical motions of the interstellar medium. No decay emission from 60 Fe 
has been observed from the Cygnus region and an upper limit of 1.6 x 10~ 5 ph cm 2 s 1 was derived. 

Conclusions. The determination in a coherent way of the gamma-ray line emission of the clustered OB population of the nearby 
Cygnus complex will allow a more accurate comparison of observations with theoretical expectations based on the most recent stellar 
models. 

Key words. Gamma rays: observations - open clusters and associations: individual: Cyg OB - Stars:early-type - ISM:bubbles - 
Nuclear reactions, nucleosynthesis, abundances 



1. Introduction 

When observing our Galaxy across the electromagnetic spec- 
trum, it is remarkable that at nearly all wavelengths, a strong 
emission arises from the Cygnus region, while most of the ra- 
diation otherwise comes from the inner Galaxy (and from a few 
close objects such as the Crab and Vela pulsars or the Cassiopeia 
A remnant). Such a strong broad-band emission points to prodi- 
gious amounts of mass and energy release in that direction, at 
intermediate distances from us. Nevertheless, the true nature, 
scale, and extent of the phenomenon underlying such a powerful 
radiative output has mostly remained unclear for a long time. 

The line of sight to the Cygnus region (assumed in the fol- 
lowing to be the part of the Galactic plane located between 
longitudes 70° and 90°) is tangential to the Local spiral arm. 
In that direction we are viewing material located close to the 
solar circle, which strongly compromises a solid evaluation of 
distances between a nd 4-5 kpc from kinematic argument^] 
(ISchneider et al.ll2006l) . Strong absorption towards the Cygnus 



1 This comes from the motions of distant objects along our line of 
sight, due to Galactic rotation, being lower than the typical proper mo- 
tions of interstellar clouds and star-forming regions (especially when 
massive star clusters are present, which is the case here) 



area (up to 32 visual magnitudes) has been another hurdle on 
the way to understan ding the layout of the region (see Fig. 1 of 
ISchneider et al.l 120061) . As a consequence, it has long been de- 
bated whether the Cygnus region results from a superposition of 
several objects/structu res distributed along the line of sight be- 
tween 0.5 and 3.5 kpc dBochkarev & SitniM1985l:IWendker et all 
1 199 It lUvaniker et alj|2001l) . or whether it is a single physical 
entity located around 1 .5 kpc and likely produced by an unusual 
concentration of stars (IVeronlll965l;ICash et al.l l"l980). The latter 
scenario was claimed to be proven by the HEAO-1 and ROSAT 
observations that revealed a large soft X-ray ring in the Cygnus 
direction, which could result from a gigantic superbubble located 
at ~ 2 kpc from us, behind the Great Cygnus Rift, and e nergised 
by a series of 30-100 supernova explosions (ICash et al.lll980h . 

iHumphrevsl (1 1978b originally inventoried 9 OB associations 
in the Cygnus region, but subsequent works have cast doubt on 
the reality of Cyg OB4, 5, and 6. At present, the clustered stel- 
lar content of the region is thought to consist of Cyg OBI, 2, 
3, 7, 8, and 9 plus a dozen open clust ers, most of which may 
well be related to Cyg OBI an d OB3 dKnodlseder et al.l 12002: 
iLe D uigou & Kn6dlsede j2002h . Most of these stellar groups are 
located between 1 and 2 kpc from us, with typical distance un- 
certainties of 0.3-0.5 kpc. Recently, a stellar census based on the 
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2MASS catalogue has revealed that the richness of the Cyg OB2 
star cluster actually exceeds the previous estimates by an order 
of magnitude. Indeed, Cyg OB2 contains about 100-120 O stars 
and could well be the first example of a young Galactic globu- 
lar cluster dKnodlsederl2000t IComeron & Pasqualil2000l) . While 
this does not constitute a definitive answer to the question of the 
origin of the strong emission from Cygnus, it seems obvious that 
such a titanic object will generate a considerable energy output 
and strong ly influences the s urrou nding interstellar medium. In 
particular, Lozinsk ava et al.l (120021) have shown that Cyg OB2, 
with its increased membership, could alone, through its collec- 
tive stellar winds, account for the X-ray superbubble and its as- 
sociated system of giant optical filaments. 

Another argument supporting the i dea of a single physica l 
entity has been provided by the work of Schneid er et alj (12006). 
From KOSMA sub-millimetre and MSX mid-infrared observa- 
tions, the authors have shown that the molecular material in 
Cygnus belongs to two huge complexes of clouds (with masses 
of 2.8 and 4.8 x 10 5 M ), which are very likely the remainder of 
a giant molecular cloud out of which Cyg OB2, OB9 and possi- 
bly OBI formed. This mass distribution is now influenced by the 
UV light and winds from these stellar associations and is a site 
of active star formation. This coherent structure lies at ~ 1.5 kpc, 
the approximate distance to the Cyg OB2 cluster. 

Although the relation of Cyg OB3, OB7, and OB8 to this 
structure has not been established, we refer in the following to 
the 6 OB associations and associated open clusters of the Cygnus 
region as the "Cygnus complex". With this term, we aim at gath- 
ering all these rich and nearby stellar groups that are relatively 
correlated in terms of distance (between 1 and 2 kpc) and that 
very likely contribute to the specificity of the Cygnus region. 
Moreover, we note that the current uncertainties on the ages and 
distances of these stellar clusters do not allow firm exclusion of 
the scenario of a connexion between them all. 

The Cygnus complex amount s to a total of ~ 170 O stars, 
among which a dozen WR stars ( Knodlse der et alj |2002D. The 
region is therefore quite young, and this is further proven by 
the lack of radio supernova remnants (SNRs) or pulsars. Indeed, 
the Green catalogu^ gives about 10 SNRs bet ween 70° and 
90°, but according to a compilation of distances by Kaplan et al. 
(2004) none of these objects can be related to an association of 
the Cygnus complex. The closest object to the Cygnus complex 
is y-Cygni (SNR 078.2+2.1), located at 1.2 kpc from us. From 
the sen sitivity of their radio survey of the region. IWendker et al.l 
(1199 11) concluded that no other SNR exists up to a distance of 
10 kpc. Additional evidence for the relative youth of the complex 
comes from the absence of pulsars tow ards Cygnus up to 4 kpc , 
as indicated by the ATNF catalogue^ (Manchester et al. 2005). 
Th ese arguments als o seem to contradict the scenario considered 
bv lCashetal.l(fl980l) in which 30-100 supernovae over the past 
few Myr would have blown the Cygnus superbubble (one should 
note, however, that the SNR ages are typically much lower than a 
Myr, which could explain the present lack; see lFrail et al 1fl994l) . 

The Cygnus region therefore constitutes a promising oppor- 
tunity to look deeper into the physics of massive stars, in par- 
ticular the nucleosynthetic processes that enrich the interstellar 
medium (ISM) with new elements. In the work presented here, 
we study the nucleosynthetic activity in the Cygnus region and 
its complex of massive stars, as traced by the gamma-ray line 



emission from the decay of 26 Al and 60 Fe. These two radio- 
isotopes are thought to be produced predominantly by massive 
stars at various stages of their evol ution, from the main seq uence 
to the supernova explosion (see Limongi & Chieffi 2006, for a 
review of the different production sites). The release of 26 Al in 
the ISM occurs via stellar winds (of Wolf-Rayet stars mainly) 
and supernova explosions, while 60 Fe is only ejected by super- 
nova explosions. With lifetimes of 1 and 2 Myr respectively, 
26 Al and 60 Fe accumulate in the ISM and their decay signal 
builds up to levels that are within reach of gamma-ray instru- 
ments. The 1809keV emission of 26 Al from the whole Galaxy 
has been imaged by CGRO/COMPTEL and INTEG RAL/SPI 
(IPluschke et alJl2001ah Ipiehl et al.ll2006l: IWang et al.ll2009h and 
shows a strong emission from the Cygnus region; the 1173 
and 1332keV Galactic emission from 60 Fe has b e en observed 
by RHESSI and INTE GRAL/SPI (Smithl l2004t lHarris et all 
120051 ; IWang et al.l l2007) but no emission from Cygnus has been 
recorded, which indicates that none or only a few supernovae 
went off and did release 60 Fe in that region. 

Because of its strong 1809keV decay emission, the Cygnus 
region was soon identified as an interesting target for study- 
ing the nucleosynthesis of massive stars. Indeed, the Galactic 
1809keV emission results from a situation involving many 
uncertain parameters: the distribution of massive stars in the 
Galactic disk and bulge, the Galactic metallicity profile, the star 
formation rate,. ..etc. Moreover, there may well be secondary 
sources of 26 Al and 60 Fe, like the supernovae of type la, the no- 
vae, or the AGB stars, with different timescales for the release 
of both isotopes. All these elements hinder accurate determina- 
tion of the expected decay fluxes to which the observations could 
be compared. In contrast, the Cygnus region and its complex of 
massive stars form a restricted, more homogeneous population, 
the characteristics of which we know better. In addition, whereas 
the Galaxy is in a steady-state situation regarding its 1809keV 
emission, the Cygnus complex is a young, evolving population, 
which thus allows assessment of not only the total yields of 26 Al 
and 60 Fe but also the history of their ejection in the ISM. 

Since the cartography of the Galactic 1809keV emis- 
sion was achieved by CGRO/COMPTEL, a number of stud- 
ies have been devoted to the nucleosynthe sis of 26 Al an d 60 Fe 



in the Cygnus region (see for examp le del Rio et al 



Pliischke et alJl2000t ICervino et al.ll2000l: IPluschke etaTl 



1996 



2002; 



Knodlsed eretal] 120021) . The latest studies agree that the ob- 



2 Green D.A., 2006, "A Catalogue of Galactic Supernova Remnants 
(2006 April version)", Astrophysics Group, Cavendish Laboratory, 
Cambridge, UK (http://www.mrao.cam.ac.uk/surveys/snrs/). 

3 http://www.atnf.csiro.au/research/pulsar/psrcat/ 



served level of the 1 809 ke V emission from Cygnus is a factor of 
2-3 above the pre dictions, when the la tter a re based on the the- 
oretica l yields by Mevnetetal. (1997) and IWooslev & Weaverl 
(fl99l and on the proven stellar content of the Cygnus re- 
gion ( including the upward revision of the Cyg OB2 member- 
ship). Pliischk e et al.l J2000) shows that the discrepancy could 
be alleviated if one includes the enhanced (but highly uncer- 
tain) yields from massive close binary s ystems in the calcula- 
tion. Alternatively, Pliis chke et alJ d2002l) argues that the prob- 
lem may well come from underestimating the stellar content 
of the Cygnus OB associations, because of a strong interstel- 
lar extinction, as was p roven to be the case for Cyg OB2. In 
iKnodlseder etal] §002), the authors eventually speculate that 
future stellar models, especially models including the effects of 
stellar rotation, could solve the dilemma. 

In recent years, the models of massive stars have in- 
deed undergone major improvements, such as the introduc- 
tion of stellar rotation or the use of revis ed mass-loss rates 
or nuclear c ross-sections (see for example iMevnet & Maeder 
2000, 12003). The effects of metal licity h ave also been ex- 
plored bv IMevnet & Maederl d2005h . while iLimongi & Chieffi 
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(2006) computed the contributions of the supernovae of type 
Ib/c in a consistent way. Meanwhile, the launch in 2002 of the 
INTEGRAL gamma-ray observatory has provided the commu- 
nity with a new instrument for addressing the problem from the 
observational side. In particular, INTEGRAL has devoted sub- 
stantial observing time to the Cygnus region through dedicated 
programmes, and its high-resolution spectrometer SPI should al- 
low the previous photometric measurements to be confirmed and 
add spectrometric information about the gamma-ray decay lines. 

The aim of our work is to re-evaluate the situation of the 
gamma-ray line emission from Cygnus, in terms of observa- 
tions, as well as theoretical predictions. For the present paper, 
we used the INTEGRAL/SPI observations to characterise in de- 
tail the 26 Al and 60 Fe decay emission from the Cygnus region. 
The analysis of the SPI data is presented in Sect.|2]and the results 
will be exposed in Sects. [3]and[4] The observational constraints 
thus obtained will then be compared in a forthcoming paper to 
theoretical predictions based on the most recent models of mas- 
sive stars' nucleosynthesis. 

2. Data analysis 

The SPI spectrometer is a coded-mask telescope with a large 
field of view of 16° x 16° and an angular resolution of 2.8° 
(detailed information about the SPI ins trument can be found in 
IVedrenne et alJ2003tlRoques et al.l2003h . And yet, the complex- 
ity of the response function means that reconstructing a given 
source intensity distribution cannot be achieved from a single 
exposure using standard correlation methods. Instead, we have 
to resort to multiple exposures (following a dithering pattern) 
and then to model-fitting approaches. 

Thanks to an array of 19 high-purity germanium detectors, 
SPI performs spectrometry of gamma radiations with energies 
between 20 keV and 8 MeV, with a spectral resolution of about 
2 keV FWHM at 1 MeV. Interactions of gamma photons with 
the detectors can be of two types: single events (SE) when the 
incoming photon deposits energy in one detector only, and mul- 
tiple events (ME) when the incoming photon deposits energy in 
several adjacent detectors through Compton diffusion and pair 
creation. In the following, we only make use of SE data. It should 
be noted that detector 2 and 17 failed on revolutions 142 and210. 

The data produced by SPI are strongly background- 
dominated (signal-to-noise ratio of less than 1 %) because sub- 
stantial hadronic interactions and electromagnetic cascades are 
induced by solar and cosmic -ray particles hitting the satellite. 
These result in an instrumental background composed of a con- 
tinuum and a certain number of lines of various intensities, 
some of which correspond to the same nuclear transitions we 
are searching for (for example the line at 1809keV from the de- 
excitation of 26 Mg inside the instrument). 

Extracting the weak gamma-ray line signals from astrophysi- 
cal objects therefore requires an accurate modelling of that back- 
ground noise; in our work, this is achieved through combinations 
of activity tracers that are hypothesised to reproduce the various 
trends and time evolutions of the background. In this effort, the 
most utilised tracer is the so-called GEDS AT, which corresponds 
to the rate of saturating events in the detectors (energy deposits 
above 8 MeV). The GEDSAT has indeed proved to be strongly 
correlated to the SPI data over a broad energy range; in addi- 
tion, its count rate (~200cts s~') is about two orders of magni- 
tude higher than the typical count rate measured at gamma-ray 
line energies, which means that it suffers little statistical uncer- 
tainty and can readily be used to draw predictions. Nevertheless, 
the GEDSAT is not a perfect tracer and it needs to be comple- 



mented and adjusted, as we see below. More information about 
the instrumental background and our methodo logy for analysing 
the SP I data can be found in lJean et alj d2003l) and lMartin et alj 
(2009), respectively. 

2.1. Data set 

The data used for this analysis are all-sky data that were col- 
lected between revolutions 7 and 494 of the INTEGRAL satellite 
(the duration of one INTEGRAL orbit is ~3 days). A filtering 
was applied to exclude all pointings during which abnormally 
high count rates (due to solar flares and periodic radiation belt 
crossings) occurred, so they preclude any detection of celestial 
signals. The total effective observation time eventually amounts 
to 63.2 Ms, among which 10.8 Ms cover the Cygnus region and 
11. 2 Ms are high-latitude pointings that we use as empty fields. 
All SE data with energies in the 1780-1830, 1150-1190, and 
1300-1 350 keV ranges were selected and binned into 1 keV-wide 
bins. 

2.2. Background models 

The physics underlying the instrumental background of the SPI 
spectrometer is complex and not fully elucidated despite con - 
siderable efforts (see for example Weidens pointner et al.l l2003). 
The development of background models based on activity trac- 
ers therefore proceeds through trial-and-error and there is no 
general prescription that works successfully over the entire en- 
ergy range of the instrument (reflecting that the origins of the 
background are specific to each energy domain). 

When scaled to the mean level of the data, these background 
models reproduce the time evolution of the instrumental back- 
ground noise reasonably well. However, as we are looking for 
gamma-ray signals that are less than 1% of the data, such a 
global scaling is not accurate enough and the background mod- 
els need to be finely adjusted. 

In the present work, we study three different gamma-ray 
lines and it turned out that each energy band required a specific 
model: 

1. 1809keV line: the background model consists of two com- 
ponents, one to account for the background continuum and 
one for the background lines. The time evolution of each 
component was assumed to be proportional to the GEDSAT 
rate in each detector, and the level of the background contin- 
uum was taken from two reference energy bands at 1786- 
1802 and 1815-1828 keV. The continuum component was 
adjusted over 1 day intervals while the line component was 
adjusted over 100 days intervals. 

2. 1 173 keV line: the background model consists of three com- 
ponents, one for the background continuum and two for the 
background lines. The time variability of the background 
continuum was assumed to follow the GEDSAT rate and the 
level of the background continuum was taken from a ref- 
erence energy band at 1 178-1 182keV. The evolution of the 
background lines intensity is given by the GEDSAT added 
to a linear increase over time (likely due to the build-up 
of medium-lived 60 Co in the instrument as a result of con- 
tinuous activation). The continuum component was adjusted 
over 1 day intervals while the line components were adjusted 
over three time periods limited by the failures of detectors 2 
and 17. 

3. 1332keV line: the background model consists of a single 
component proportional to the GEDSAT rate and adjusted 
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over 3 days intervals. In this energy range, the variabil- 
ity of the background noise was stronger and thus required 
to adjust the model on shorter timescales, which decreased 
the sensitivity of the instrument to the weak celestial sig- 
nal. As a consequence, the 1332 keV line measurement only 
marginally contributes to the 60 Fe results. 



The background models presented above are fitted to the data, 
with the proper parameter set, simultaneously to models for the 
celestial contribution. 



2.3. Model fitting 

As already mentioned, the SPI data cannot be inverted and the 
astrophysical signals are extracted by fitting models of intensity 
distributions to the data. This can be done with iterative methods 
such as Richardson-Lucy or Maximum Entropy that reconstruct 
the sky intensity distribution pixel by pixel, by setting the flux in 
each sky element to the value that best accounts for the data (in a 
statistical sense). An alternative is to set the shape of the source 
model (from observations at other wavelengths for instance) and 
to adjust its overall flux. The latter method is the one we used be- 
cause it ensures that the result will be astrophysically meaning- 
ful and compatible with the performances of the instrument (on 
angular resolution especially), while reconstruction algorithms 
often generate artefacts from statistical fluctuations and uneven 
exposure. 

In the present work, simple analytical models, such as two- 
dimensional Gaussian profiles and a ring model for the Galaxy, 
were used for the sky contributions (the details of these mod- 
els will be given below). Each of these intensity distributions is 
mapped into the data space through a convolution with the in- 
strument response function (IRE) for each pointing of the data 
set and then forms a sky model component. The sky model 
components and the background models are then fitted simul- 
taneously to the data using a Maximum Likelihood criterion for 
Poissonian statistics. 

Once all models (sky and background) were fitted to the 
data, the quality of the analysis was assessed by examining the 
residuals (data counts minus fitted background and sky models 
counts). Since we are working with background-dominated data, 
the modelling of the instrumental background can easily intro- 
duce systematic errors that might be compensated for by artifi- 
cial signals from the sky. In addition, the assumption made on 
the sky intensity distribution might well be inappropriate or in- 
complete, which would introduce bias or leave unaccounted for 
events in the data space. 

To evaluate the quality of our analysis, we studied the 
distribution of the residuals after back-projection on the sky. 
Basically, for all pointings of the data set, the residuals in each 
detector were uniformly projected across the mask (or more pre- 
cisely "across" the instrument response function) over the sky re- 
gion intersected by the field of view. The distribution of these sky 
residuals should then follow a statistical distribution, the shape 
of which is derived from simulat ed observa t ions. T his method 
is described in greater detail in iMartin et al.1 (120091) and proved 
that our hypotheses regarding both the background noise and the 
sky intensity distribution are satisfactory and that the results ob- 
tained are not affected by systematic errors. 



3. The gamma-ray line emission from Cygnus 

3.1. Morphology of the 1 809 keV emission 

We have characterised the morphology of the 1809keV emis- 
sion from the Cygnus region through a model-fitting approach. 
Different models, in terms of size and position over the Cygnus 
region, were fitted to the data (simultaneously to a model for the 
instrumental background), and after a correction for the num- 
ber of independent trials, the model with the highest likelihood 
criterion was selected. This method has been tested on simu- 
lated observations with various true source sizes and positions 
and proved successful. 

Due to small statistics on the astrophysical fluxes and poor 
angular resolution of the instrument, the morphology of the 
Cygnus 1809keV emission cannot be determined with a high 
degree of accuracy. Therefore, we focused on determining of the 
extent and centre of the 1809keV emission. We considered 6 
model types: point-source and 2D Gaussians with standard de- 
viations cr of 1°, 2°, 3°, 4°, and 5° (the typical size of these 
Gaussian models is given by ~ 3<x). Each of these 6 models 
was moved over a grid of 176 positions ranging from 65° to 95° 
in longitude and -10° to 10° in latitude, with 2° steps in each di- 
rection. 

From this analysis, the best-fitting model is a 3°x3° Gaussian 
centred on (l,b) = (81°, -1°). Typical uncertainties are ~ 3° on 
the size of the model and ~ 3° on the position (in both longitude 
and latitude). However, the maximal extent of the diffuse emis- 
sion at 1809keV from Cygnus is poorly constrained because of 
intrinsic limitations of the coded-mask imaging system. Indeed, 
a coded-mask instrument is unable to detect a flat intensity distri- 
bution because the latter leads to a uniform flux over the detector 
plane, which cannot be distinguished from instrumental back- 
ground. Flat or very extended intensity distributions can only be 
detected by alternating pointings on the source (ON pointings) 
with pointings to empty fields (OFF pointings), in which case 
the coded-mask instrument is used as a mere collimator. For the 
Cygnus extended emission at 1 809 ke V, it turned out that sources 
with typical sizes above 8-10° cannot be distinguished (provided 
most of the emission remains within the field of view); therefore, 
our best-fitting model should formally be considered as a lower 
limit. It should be noted, however, that the typical size of the 
Cygnus 1809 keV emission on the COMPTEL map is consistent 
with our best-fitting model. 

3.2. Characteristics of the 1 809 keV emission 

By representing the 1809keV emission from the Cygnus region 
by the 3° x 3° Gaussian model determined in l3.ll the fit to the 
data at 1780-1830keV yields the spectrum of Fig. Q] The cen- 
troid of the line is 1808.8 + 0.4 keV, which is consistent with 
the laboratory value of 1808. 6keV and hence indicates that the 
emitting medium has no or only a small bulk motion with respect 
to the observer. The width of the line is 3.4 + l.OkeV (FWHM) 
and is consistent with the 3.0 keV (FWHM) spectral resolution 
of the SPI spectrometer. Considering the statistical uncertainty 
on the line width, the emitting medium may be affected by ex- 
pansion or turbulent velocities with values up to ~ 200kms~ 1 . 

Based on our Gaussian model, the flux from the Cygnus re- 
gion in the 1 806- 1 8 1 2 keV band is (6.0 ± 1 .0) x 1 0~ 5 ph cm 2 s" 1 . 
This value, however, gives the total flux seen towards the Cygnus 
region and thus includes the contribution from the Cygnus com- 
plex located at about 1.0-2.0 kpc from the Sun, plus foreground 
and background contributions from the Galaxy as a whole. If a 
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Fig. 1. Spectrum of the 1809keV emission from the Cygnus re- 
gion, from about 4 years of INTEGRAL/SPI observations. The 
red line represents the best Gaussian fit to the data points. 



Fig. 2. Spectrum of the combined 1 173/1332keV emission from 
the Cygnus region, from about 4 years of INTEGRAL/SPI ob- 
servations. The red line represents the best Gaussian fit to the 
data points. 



model for the Galactic 1809keV emission (to be described in 
14. U is fitted to the all-sky data simultaneously to our Gaussian 
model, the flux associated with the Gaussian model falls to 
(3.9 + l.l)xl0~ 5 phcrrT 2 s _1 . We interpret the latter flux as being 
mostly caused by the OB associations and open clusters of the 
Cygnus complex, while the remaining ~ 2.0 x 1CT 5 phcirr 2 s _1 
are very likely due to a more diffuse population, distributed all 
along the line of sight. This population may include, for instance, 
unclustered WR and OB stars (runa way stars or B stars from pre- 
vious star formation episo des, see IComeron & Pasquaul d2007l) 
and lComeron et al.l (12008)). or isolated SNRs. 

3.3. Search for 1 1 73/1 332 keV emission 

As is shown in I4.21 the Galactic 60 Fe decay emission is about 
7 times weaker than the Galactic 26 Al decay emission. The 
1 809 ke V line allows morphology studies, either over the Galaxy 
or more specifically in the Cygnus region, which cannot be en- 
visaged with the 60 Fe lines. In the case of the Cygnus region, a 
study of the morphology of the 1 173/1332keV emission may be 
impossible simply because there may not be any 1 173/1332keV 
signal at all. Indeed, since the region seems quite young, it may 
well be that none or few supernova occurred and thus little 6l) Fe 
has been released in the ISM, if any. 

Using the same model for the Cygnus region as in 13.21 and 
fitting this model to the data at 1150-1190 and 1300-1350keV 
gives the spectrum of Fig. [2] No signal shows up and an upper 
limit on the flux of 1.6 x 10~ 5 phcirr 2 s~' can be derived (2<x 
upper limit assuming the line is not broadened, as is the case at 
1809keV). We emphasise here that this upper limit applies to 
the Cygnus region and is therefore conservative if used for the 
Cygnus complex. 

3.4. Discussion 

The characteristics of the gamma-ray line emission from Cygnus 
agree with the current picture of the region as derived from ob- 
servations at other wavelengths. 

First of all, the absence of 60 Fe decay signal from that di- 
rection suggests that no or very few supernovae went off in the 
Cygnus complex, which is compatible with the lack of SNRs or 



radio pulsars. Then, the extended emission from Al decay is 
centred on a position which is consistent with the position of the 
Cyg OB2 cluster; the latter is expected to dominate the energet- 
ics and nucleosynthesis of the complex given its richness com- 
pared to the other OB associations. The negligible shift in the 
1809 keV line is in accordance with the expected relative motion 
due to Galactic rotation for the direction and mean distance of 
the Cygnus complex, and the width of the line is consistent with 
typical ISM motions. 

The extent of the 1809keV emission reflects the diffusion of 
the ejected 26 Al in the ISM; in particular, the 26 Al is expected 
to fill to some extent the superbubbles blown by the various OB 
associations. Inside these low-density cavities, the radioisotope 
is affected by expansion motions that accompany the growth of 
the superbubble and may cause some broadening of the decay 
line. 

Our estimate of the total 1809keV flux from the Cygnus re- 
gion is (6.0 ± 1 .0) x 10~ 5 ph crrr 2 s , in agreement with previous 
measurements. The flux values obtained with COMPTEL were 
in the range 7-8 x 10~ 5 ph crrT 2 s , depending on the integration 
region and imaging method used, with a ty pical uncertainty of 
about 1.5 x!0~ 5 phcm ' 2 s _1 (see for instance ldel Rio et alJ l996; 
iPliischke et al]l20~01bl) . Other measurements from SPI are in the 
range 6-7 xl0~ 5 phcirr 2 s , here again depending on the imag- 
ing method and integration r egion, with a typical uncertainty of 
about 1.2 xl0~ 5 phcirr 2 s" 1 (IWangll2007l) . 

The most interesting result of our study is that analysing 
the Cygnus region within a Galactic context has revealed that 
the 1809keV flux really attributable to the Cygnus complex of 
massive stars is only ~ 65% of the total flux coming from that 
direction. This point is essential to the comparison of observed 
fluxes with the theoretical predictions, because the latter are of- 
ten mostly based on the clustered stellar content, which is better 
inventoried than the isol ated or non-clustered sources. It is inter- 
esting to note here that IPliischke et all d2002l) find that the con- 
tribution of non-clustered sources may amount to about 20-30%, 
which tallies with our observational estimate; however, their re- 
sult was based on an assumed correction of the stellar content of 
the Cygnus OB associations that has never been observationally 
substantiated afterward. 
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Last of all, we want to emphasise that our result is based 
on the assumption of a Gaussian profile for the emission of the 
Cygnus complex, whereas the true intensity distribution may 
well be more structured. As already said, the SPI (or even 
COMPTEL) data do not allow accurate determination of the 
morphology of the diffuse 1809keV emission from Cygnus and, 
as such, our assumption of a Gaussian profile for the extended 
source is a potential source of error. 

4. The gamma-ray line emission from the Galaxy 

4.1. The 1 809 keV Galactic emission 

The detailed study of the Cygnus gamma-ray line emission re- 
quired to properly take the mean Galactic emission into account 
in that direction. Therefore, the following presents the model 
used for the mean Galactic emission at 1809keV. 

When studying the 26 Al Galactic emission, the most impor- 
tant quantities to be determined are the total mass of 26 Al and 
its distribution throughout the Galaxy; and yet, as already stated, 
the small statistics on the astrophysical fluxes and the poor an- 
gular resolution of the instrument severely restrict the informa- 
tion that can be derived from the SPI observations. The main 
source of Galactic 26 Al are the massive stars, the vast majority 
of which clusters in the spiral arms, but our data do not have the 
degree of accuracy necessary for unambiguous reconstruction 
of a sophisticated spiral arm pattern from the intensity distribu- 
tion observed on Earth. Instead, we focus on global quantities 
such as the azimuthally-averaged radial distribution of 26 Al in 
the Galaxy. 

We used a five-component model in which the Galaxy is di- 
vided into five concentric rings, each ring containing a certain 
mass of 26 Al distributed as 

P (Ri < R < R M ) = Ki x e'i (i=l,...,5) (1) 

where R is the galactocentric radius, the values define the 
edges of the rings, and the Kj are the scaling factors set by the 
model-fitting process. The distribution of the radio-isotope is as- 
sumed to be uniform radially inside each ring and exponential 
vertically, with a characteristic height zo- We started our analysis 
with a value of zq = 90 pc, which corresponds to the character- 
istic thickness (within the solar circle) of the molecula r gas disk 
out of which massive stars form ( Bronf man et al.l 2000). This pa- 
rameter will be discussed in more detail later on. 

This five-component model turned out to be insufficient to 
account for the data and should be supplemented by two com- 
ponents: the Gaussian model determined in !3. ll for the Cygnus 
diffuse emission and another Gaussian model for a very extended 

Table 1. Galactic mass distribution of 26 Al obtained from our 
analysis of 4 years of INTEGRAL/SPI observations (see text). 



ring 1 (0-3 kpc) 


0.17 ±0.07 Mo 


ring 2 (3-6 kpc) 


1.18 +0.15MQ 


ring 3 (6-8 kpc) 


0.18 ±0.08 Mq 


ring 4 (8-9 kpc) 


0.15 ±0.04 Md 


ring 5 (9-12 kpc) 


-0.28 ±0.18M o 


Total 26 Al mass 


1.68 ±0.19Md 



emission from the Sco-Cen region. The Sco-Cen emission was 
found to be represented well by a 5° x 5° Gaussian model cen- 
tred on (I, b) = (0°,30°), which is nearly consistent in position 
and extent with the superbubble b lown by the Upper-Scorpius 
OB association (see lde Geiisll 19921) . 

Fitting these seven intensity distributions (five Galactic rings 
+ Cygnus + Sco-Cen) to the data (simultaneously to a model 
for the instrumental background) gives the 26 Al mass distribu- 
tion summarised in Table [TJ (Note that ring 5 is not used in the 
computation of the total 26 Al mass because of its non-physical 
value and large uncertainty). The total corresponding flux from 
the Galactic plane in the 1806-1812keV band amounts to (8.4 + 
1 .7) x 10~ 4 ph cirT 2 s _1 . In addition to this azimuthally-averaged 
radial distribution of 26 Al in the Galaxy, the nearby regions of 
Cygnus and Sco-Cen bring additional emission at a level of (3.9 
+ 1.1) and (6.2 + 1.6) x 10~ 5 phcm~ 2 s _1 respectively. This pa- 
per focuses on the Cygnus region and the 1809keV e mission 
from Sco-Cen will be discussed in a forthcoming paper (Diehl 
120091) . 

The spectrum of the 1809keV emission from the Galactic 
disk (excluding the Sco-Cen and Cygnus contributions) is shown 
in Fig. [3] The centroid and width of the line are consistent with 
no Doppler shift and/or broadening, but it should be emphasised 
that this spectrum results from the contributions of many re- 
gions along the Galactic plane, each region having specific phys- 
ical properties like the proper motion (relative to us) due to the 
Galactic rotation pattern. A study of the spectral characteristics 
of the 1809 keV l ine along the Galactic plane can be found in 
iDiehl et all l2006h and lWang et alJ d200l . 



4.2. The 11 73/1 332 keV Galactic emission 

Due to lower intensity, the decay emission of 60 Fe cannot be 
analysed in as much detail as the emission from 26 Al; in par- 
ticular, the spatial distribution of the 60 Fe isotope in the Galaxy 
cannot be estimated, even for a simple model. The only observa- 
tional constraint accessible to us is the emission flux, especially 
the flux from the Galactic centre region (where most of the emis- 
sion and exposure time are concentrated). 

The decay emission at 1 173 and 1332 keV from the 60 Fe 
isotope is assumed to follow the distribution of massive stars 
throughout the Galaxy (it should be remembered here that, in 
contrast to 26 Al, it has never been demonstrated observation- 
ally that massive stars are the dominant source of 60 Fe). As a 
tracer of the massive stars, we used the infrared emission from 
the interstellar dust heated by the strong UV luminosity of OB 
stars, as observed at 240 /mi by the COBE/DIRBE instrument. 
Fitting this intensity distribution to the data in the 1 1 50-1 190 and 
1 300-1 350 keV bands (simultaneously to a background model 
for each energy range) and combining the resulting spectra gives 
the spectrum of Fig.[4]for the Galactic 60 Fe decay emission. 

The line is too weak to allow for a solid spectrometric study, 
so we focus on photometric considerations. Assuming the line 
width is purely instrumental, the flux in the 1 173 and 1332keV 
decay lines is (9.1 + 2.6) x 10~ 5 phcrrT 2 s _1 for the whole 
Galaxy and (3.5 + 1.0) x 10~ 5 phcirT 2 s _1 for the central re- 
gion (defined as a steradian centred on the Galactic centre). It is 
worth mentioning that the flux value obtained for the Galactic 
centre region is very likely more reliable than the one for the 
whole Galaxy. Indeed, due to the weak emission intensity in the 
Galactic disk and because of the strongly inhomogeneous expo- 
sure of INTEGRAL, the scaling of the DIRBE intensity distribu- 
tion is driven by the Galactic centre region; as a result, the flux at 
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high longitudes is more a consequence of this "biased" scaling 
than a solid measurement. 



4.3. Discussion 

The chacteristics of the 26 Al and 60 Fe decay emissions from 
the Galaxy agree with the results obtained in previous stud- 
ies using INTEGRAL/SPI or other instruments. The flux values 
are consistent with all previous measurements. Our analysis of 
the 1173/1332keV lines yielded a flux in the inner Galaxy of 



(3.5 
(3.6 



+ 1.0) x 10" 
± 1.4), (3.7 : 



5 ph cm 2 s 1 , in agreement with the values of 
: 1.1 ) and (4.4 + 0.9) x lO ^phcm^s" 1 pro- 
duced from RHE SSI (lSmithll2004l) and SPI (lHarris et al.ll2005l 
IWang et aTll2007l) observations. 

Regarding the Galactic 1809keV emission, the emission in- 
tensity we d erived is also consistent with previous estimates (see 
for instance iDiehl et all 120061 IWang et all [20091) . We found a 
26 Al stationary mass of 1.7 +0.2 M (excluding the specific Sco- 
Cen and Cygnus contributions), most of which is confined be- 
tween 3 and 6kpc, at the position of the Giant Molecular Ring. 
This value, however, is linked to an assumption on the latitudinal 
extent of the 1809keV Galactic emission. This parameter cannot 
be strongly constrained by the SPI observations due to intrinsic 
limitations of the coded-mask imaging system (see the detailed 
explanation given in !3. lb . Therefore, we adopted a characteristic 
value of 90 pc for the vertical extent of the 1809keV emission. 
This is a "minimal" hypothesis since it is the characteristic scale 
height of the disk of molecular gas out of which massive stars 
form; since massive stars are the dominant source of 26 Al, the 
vertical extent of the 1809keV emission cannot be below these 
90 pc. Nevertheless, the ejected 26 Al will diffuse in the ISM, in 
particular inside the superbubbles blown by the collective winds 
and SN explosions of the OB associations, and can travel out to 
distances of a few hundred pc. The 1809keV emission is there- 
fore expected to be more extended vertically than 90 pc. To eval- 
uate the importance of this parameter, we performed the analysis 
of the 1809keV all-sky emission using as a disk height for the 
Galactic ring model a value of 300 pc instead of the minimal 
90 pc. The fit to the data (using the same background model) 
was found to be identical in terms of statistical criterion, but the 
Galactic mass of 26 Alroseto2.0±0.2M Q . 

Our estimate of the 26 Al stationary mass in the Galaxy is 
lower than but compatible with, w ithi n the error bars, the values 
put forward by Die hl et all (120061) and lWang et all (120091) . which 
are respectively 2.8 +0.8 M and 2.7 ±0.7 M . The larger un- 
certainties associated with these values actually reflect our poor 
knowledge of the exact sky distribution of the 1809keV emis- 
sion and associated three-dimensional Galactic distribution of 
the 26 Al isotope. Indeed, these authors have used, as we did in the 
present work, a model-fitting approach to extract the photomet- 
ric (and spectrometric) information about the 1809 keV Galactic 
emission; in order to evaluate the dependence on the assumed 
galactic distribution, they have considered a series of plausible 
models, and the above-cited values reflect the scatter in the in- 
ferred 26 Al masses. Our work complements this study in that our 
concen tric ring model has not been considered by Dieh l et aT] 
(2006|) and IWang et all d2009lh and this model implies a 26 Al 
mass at the lower end of the mass range they determined. On 
th at point, it is i ntere sting to note that, among the models tested 
by I Wang et all (120091) . the best fitting ones are actually associ- 
ated with the lowest 26 Al mass values, with a meaningful (but 
not drastic) statistical improvement (Michael Lang 2008, private 
communication). A major difference in our analysis is also that 
we separately modelled the Cygnus and Sco-Cen complexes of 
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Fig. 3. Spectrum of the 1809keV emission from the Galaxy 
(Sco-Cen and Cygnus regions excluded), from about 4 years of 
INTEGRAL/SPI observations. The red line represents the best 
Gaussian fit to the data points. 
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Fig. 4. Spectrum of the combined 1 173/1332keV emission from 
the Galaxy, from about 4 years of INTEGRAL/SPI observations. 
The red line represents the best Gaussian fit to the data points. 



massive stars. It may be that these peculiar nearby regions, on 
which there is considerable SPI exposure, somehow drive the 
fit of some all-sky models to higher underlying 26 Al masses. 
Ongoing INTEGRAL observations of the Galaxy, and specifi- 
cally of the Cygnus and Sco-Cen regions, will help refining our 
perception of the 1809keV emission and associated 26 Al distri- 
bution. 

As already said in 13.21 a major result of our all-sky analy- 
sis of the 1809keV signal is a better estimate of the 1809keV 
decay flux attributable to the Cygnus complex. We indeed disen- 
tangled the contribution of the nearby clustered OB population 
of the Cygnus complex in a coherent way from the contribution 
of the more diffuse and spatially extended non-clustered popula- 
tion likely to contribute to the emission. Also noteworthy is the 
non-negligible 1809keV flux found from the Sco-Cen region. A 
similar flux was found from a lmost the same data set and is fur- 
ther discussed bv lDiehll d2009h . 
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5. Conclusion 

We have analysed about 4 years of data from the 
INTEGRAL/SPI spectrometer to characterise the 26 Al and 
60 Fe decay emission from the Cygnus region in detail, with 
the purpose of constraining the nucleosynthesis activity of the 
Cygnus complex of massive stars located between 1.0 and 
2.0 kpc from the Sun. 

The 26 Al decay emission from Cygnus can be represented 
as a 3° x 3° Gaussian centred on (l,b) = (81°, -1°), a position 
consistent with that of the massive Cyg OB2 cluster thought 
to dominate the energetics and nucleosynthesis of the Cygnus 
complex. The total 1809keV flux from the Cygnus region is 
(6.0 + 1.0) x 10~ 5 phcirT 2 s~' and the flux really attributable 
to the Cygnus complex, obtained by properly considering the 
mean Galactic background and foreground contribution, reduces 
to (3.9 + 1.1) x 10~ 5 phcirT 2 s _1 . This effort to coherently disen- 
tangle the contribution of the rather well-identified clustered OB 
population of the nearby Cygnus complex from the contribution 
of the more diffuse and less known non-clustered population 
will allow a more accurate comparison of observations with 
theoretical expectations based on the most recent stellar models 
(to be presented in a forthcoming paper). 

The 1809keV line centroid agrees with expectations con- 
sidering the direction and distance of the Cygnus complex, and 
the line width is consistent with typical ISM motions. No decay 
emission from 60 Fe was observed from the Cygnus region and 
an upper limit of 1.6 x 10~ 5 phcirT 2 s _1 was derived. 

In the course of our work, a modelling of the Galactic 26 Al 
content and decay emission was performed. In agreement with 
previous results, we have found a total 1809keV flux of (8.4 + 
1.7) x 10~ 4 phcirT 2 s -1 . Our model implies a 26 Al stationary 
mass of 1.7-2.0 ±0.2 M Q , depending on the ill-constrained 
vertical extent of the 1809keV Galactic emission, and most of 
this mass appears to be located in the Giant Molecular Ring, 
between 3-6 kpc. We also observed a weaker Galactic 60 Fe de- 
cay emission with a total flux of (9.1 + 2.6) x 10 _5 phcm _2 s _1 , 
which reduces to (3.5 + 1.0) x 10~ 5 phcirT 2 s~' in the central 
region of the Galaxy. 
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